Drift-Controlled Anomalous Diffusion: A Solvable Gaussian Model 
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We introduce a Langevin equation characterized by a time dependent drift. By assuming a temporal 
power-law dependence of the drift we show that a great variety of behavior is observed in the 
dynamics of the variance of the process. In particular diffusive, subdiffusive, superdiffusive and 
stretched exponentially diffusive processes are described by this model for specific values of the 
two control parameters. The model is also investigated in the presence of an external harmonic 
potential. We prove that the relaxation to the stationary solution is power-law in time with an 
exponent controlled by one of model parameters. 
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Diffusive stochastic processes, i.e. stochastic processes 
x{t) characterized by a linear grows in time of the vari- 
ance < x'^{t) >cx t are quite common in physical sys- 
tems. However deviations from a diffusive process are 

\ observed in several stochastic systems. Superdiffusive (< 
x^{t) >(x with V > I) and subdiffusive (< x'^{t) >cx 
with V < 1) random processes have been detected and in- 
vestigated in physical and complex systems. A classical 
example of superdiffusive random process is Richardson's 
observation that two particles moving in a turbulent fluid 
which at time i = are originally placed very close the 
one with the other have a relative separation i at time t 
that follows the relation < £'^{t) >cx Most recent 

[ examples include anomalous kinetics in chaotic dynamics 
due to flights and trapping [^[jsj, anomalous diffusion in 
aggregate of amphiphilic molecules and anomalous dif- 
fusion in a two-dimensional rotating flow Q . Subdiffusive 
stochastic processes have also been detected and investi- 
gated. Examples includes charge transport in amorphous 
semiconductors [^,^ and the dynamics of a bead in poly- 

I mers . Another class of stochastic processes which are 
not diffusive in a simple way is the one characterized by a 
variance with a stretched exponential time dependence. 
When a such process is Gaussian distributed the prob- 
ability of return to the origin Po{t) is described by the 
Kohlrausch law Po{t) oc exp[— i"^] with i/ < 1. Similar 
behaviors are observed in glassy systems and in random 
walks in ultrametric spaces ^. 

The modeling of some of the above discussed anoma- 
lous diffusing stochastic processes has been done by using 
a variety of approaches. To cite some examples, we recall 

' that superdiffusive and subdiffusive processes have been 
modeled by writing down a generalized diffusion equation 
|^, P^|Jfl}| , by introducing Levy walks models |fj] , by using 
a fractional Fokker-Planck equation approach |l3 or by 
using "ad hoc" stochastic models such as, for example, 
the fractional Brownian motion p^ . 

In this rapid communication we introduce a class of 
Langevin equations able to describe all the different 



anomalous regimes discussed above for Gaussian pro- 
cesses. Specifically, we study the properties of the class 
of Langevin equations 

x + j{t)x^T{t), (I) 

where j{t) is a function of time t and T(t) is a Langevin 
force with zero mean and with a correlation function 
given by < T{t2)T{ti) >= DS{t2 - h). Equation (f) 
describes an Ornstein-Uhlenbeck process []f5| in the par- 
ticular case of 7(t) — 7. This equation is linear and 
solvable. For the sake of simplicity we set the boundary 
condition of Eq. (1) at t = 0. The formal solution of Eq. 
(1) is 

x{t)^xiO)G{t) + G{t) f^.ds, (2) 
Jo G(s) 

where G(t) = exp[— j{s)ds]. By using this formal so- 
lution and all order correlation functions of T(t) we ob- 
tain all central moments of x(t). The first two central 
moments are given by 

< x{t) >= x{0)Git) = ^{t), 

< ix{t) - ^,{t)f DG^it) EE o\t). (3) 

The general relation between higher-order even central 
moments and the second central moment of the investi- 
gated processes is the one observed in a Gaussian pro- 
cess. Moreover odd central moments are zero hence we 
conclude that the stochastic processes described by Eq. 
(f) is Gaussian. 

We now consider the two-time correlation functions of 
the process x{t) and of its time derivative x{t). In the 
following we label the two times ti and t2 of the corre- 
lation functions in such a way that t2 > ti. By using 
the formal solution of Eq. (2) we determine the two-time 
correlation function for the random variable x{t) 



< X{ti)x{t2) >= H{h)n{t2) + 



G{t2 



(4) 



In general, the correlation function < x{ti)x(t2) > is not 
a function of t2 — ti and therefore the process is usually 
non-stationary. 

By starting from the correlation function of x{t) and 
from the formal solution of the Langevin equation we 
obtain the two-time correlation function of x{t) as 

< X{ti)x{t2) fiv{tl)f^v{t2) + Dd{t2 - h) 

+l{t2)G{t2)F{h), (5) 

where ^.^{t) = —^{t)fi{t) indicates the mean of the time 
derivative x{t) and F{ti) = [y{ti)a'^{ti) - D] /G{ti). 
The two-time correlation function of x{t) is the sum of a 
delta function and of a smooth function. 

The Fokker-Planck equation associated to the 
Langevin equation given in Eq. (1) is 



dp d D d^p 

- = -(7(tM + -^. 



(6) 



This Fokker-Planck equation is the same as the Smolu- 
chowski equation of a Brownian particle moving in a har- 
monic oscillator with a time dependent potential U (x) cx 
x'^j{t). In our study we consider both positive and nega- 
tive values of j{t). For positive values of 7(f) the position 
X = is a stable equilibrium position whereas in the op- 
posite case X = is an unstable equilibrium position. 

We calculate the two-time conditional probability den- 
sity P{x2,t2\xi,ti) as the Green function of the Fokker- 
Planck equation. In our formalism ti < t2. Our determi- 
nation is done by working with the Fourier transform of 
P{x2,t2\xi,ti) with respect to the x variable. The equa- 
tion for the Fourier transform of P{x2,t2\xi,ti) is a first 
order partial differential equation, which can be solved 
by the methods of characteristics. We obtain 



P{x2,t2\xi,ti 



1 



exp 



where m{t2,ti) 



y/2Trs^{t2,h) 

{X2 " 771(^2, 

2s2(i2,ti) 



(7) 



It^ liy)dy , and s'^{t2,ti) = 



exp 

abihty of Eq. (7) is 
bility. Moreover, Eq. 
Kolmogorov equation, 
tegration one can verify that 

/ P{x3,t3\x2,t2)Pix2,t2\xi,ti)dX2. 

In the rest of this rapid communication we restrict our 
attention to the class of Langevin equations with a drift 
term which has a temporal behavior of the form 



dz. Hence the transition prob- 

a Gaussian transition proba- 
(7) satisfies the Chapman- 
In fact from a direct in- 

P{x3,t3\xi,ti) = 



7(0 ~ a/t^ 



(8) 



for large time values. We study the stochastic process 
of Eq. (1) for different values of parameters a and /3. 



Specifically we focus on the asymptotic temporal evo- 
lution of the variance and of the two-time correlation 
function of x{t). We recall that for a — Eq. (1) de- 
scribes a Wiener process with a variance increasing in 
a diffusive way, (T^(<) ~ t, and a delta correlated x{t). 
When /3 = Eq. (1) describes an Ornstein-Uhlenbeck 
process and two regimes are observed depending on the 
sign of a. When a > the stochastic process has a sta- 
tionary Gaussian solution, whereas when a < there is 
no stationary state and the variance increases asymptot- 
ically in an exponential way, a{t) ~ exp(2|a|i) p^ , p^ . 
The two-time correlation of the velocity decreases in an 
exponential way as exp(— |a|(t2 — ti)). 

The cases considered above are known, in addition to 
these cases we observe a large variety of new behaviors 
controlled by the specific values of parameters a and /3. 
By investigating the (/3, a) set of parameters, we detect 
different anomalous behavior that we discuss below sys- 
tematically by considering different regions of the /3 pa- 
rameter. 

(i) Region with /3 > 1. The process x{t) is diffusive and 
its variance increases linearly with time for any value of a. 
The two-time correlation function of can be obtained 
starting from Eq. (5). A direct calculation gives 



< x{ti)x{t2) >- 



■ exp 



—at 



1-/3 



F{ti). 



(9) 



The process x{t) can be positively or negatively corre- 
lated depending on the sign of a. When a > (a < 0) 
the correlation is negative (positive). This property is 
valid for any value of P parameter. By investigating the 
explicit form of Eq. (9) one observes that the correlation 
function decreases as a function of t2 with a power-law 
dependence, < x{ti)x{t2) >~ l/if- 

(ii) Region with P — 1. In this case we observe two 
regimes. When a > —1/2 the variance increases in a 
diffusive way cr^(t) ~ t. We find a different behavior 
when a < —1/2. In fact, by using Eq. (3) one can show 
that 



2\a\ 



(10) 



Therefore the particle performs a Gaussian superdiffusive 
random process. At the boundary value a = —1/2 the 
variance increases in a log divergent way as {t) ~ t log t. 
The two-time correlation function of x(t) is determined 
starting from Eq. (5). An explicit calculation gives 



< x{ti)x(t2) 



,1+. 



(11) 



The two-time correlation function of x{t) shows a power- 
law time dependence and the x{t) process is a strongly 
dependent random process ||l^]. We wish to point out 
that when a = — 1 the diffusion of the x{t) process is 
ballistic. This specific case has been already investigated 
by E. Nelson in the framework of stochastic mechanics. 
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The Ito equation describing the stochastic process associ- 
ated with the free evolution of a Gaussian quantum wave 
packet is p8| 



t — c 

dx{t) — ^5 7tx dt + dw{t), 

-f- 



(12) 



where w{t) is a Wiener process and c is a constant . This 
stochastic equation describes the same random process 
of Eq. (1) for large values of t. 

(iii) In the region < /? < 1 we observe two regimes, 
which depend on the sign of a. When a > the variance 
increases as 



(13) 



This behavior is the customary behavior observed in sub- 
diffusive random process. The two-time correlation func- 
tion of x{t) behaves asymptotically as 



< x{ti)x{t2) >- 



aD 



■ exp 



1-/3 



(14) 



If the time interval t = t2 — ti is shorter than ti 
the power-law term dominates in this equation and 
the stochastic process is power-law anti-correlated. For 
T >> ti the two-time correlation function decreases ex- 
ponentially. 

When a < the variance increases as a stretched ex- 
ponential 



a^{t) ~ exp 



2|q| , 

1-/3 



1-/3 



(15) 



Since the process is Gaussian, the probability of return 
to the origin follows the Kohlrausch law, p{x{0),t) = 
l/^/2^Ta{t) ~ exp [— '^]- This kind of anomalous dif- 
fusion has been observed in glasses and in random walks 
on an ultrametric space The two-time correlation 
function of x{t) increases with time as Eq. (14). 

(iv) Region /? < 0. This region is essentially different 
from the previous ones because the absolute value of the 
drift term increases in time and eventually diverges. In 
this case we also observe two regimes depending on the 
sign of a. When a < the time evolution of the vari- 
ance is formally the same as Eq. (15) of case (iii). In 
this region of P parameter the variance increases more 
than exponentially in time. When a > we find that 
the variance decreases with time with a power-law de- 
pendence (i) ^ 1/^1/31, By using the Smoluchowski 
picture, we can interpret this behavior as the motion of 
Brownian particle moving in a time dependent potential 
which leads to a localization of the particle in the point 
X = 0. For both regimes the two-time correlation func- 
tion of x{t) is given by Eq. (9). 

We summarize the above discussed variety of diffusive 
behavior of cr'^{t) in Table I. 



The results obtained above refer to x{t) random pro- 
cesses which are not stationary. We now consider the 
problem of a process x{t) whose dynamics is controlled 
by a modified version of Eq. (1) in which the effects of 
the presence of an "external" time-independent poten- 
tial are taken into account. To this end we consider the 
specific case of a overdamped particle of mass M mov- 
ing in a viscous medium in the presence of a potential 
having a time-dependence of the kind of Eq. (8) and a 
time- independent part. The equation of motion of such 
a system is 



Mrii + g{t)x - F{x) = MT{t), 



(16) 



where rj is the friction constant and T{t) is the Langevin 
force with diffusion constant 2'qkBT /M . This equation 
is formally equivalent to 



x + -i(t)x + V\x) ^T(t), 



(17) 



when V'{x) = -F{x)/Mf], j{t) = g{t)/Mr] and D = 
2kBT /Mrj. The prime in V{x) indicates spatial deriva- 
tive. It is worth pointing out that when ^(t) goes to zero 
as t increases (as, for example, in the case "f{t) ~ a/t^ 
with > 0) Eq. (17) might have a stationary solution. 
The presence of a stationary solution depends on the ex- 
act shape of V{x). 

To investigate in a concrete example the relaxation dy- 
namics of the probability density function of x{t) towards 
the stationary solution we study Eq. (17) in the pres- 
ence of an external harmonic potential, V{x) = ^kx^. In 
this case the process has a stationary state. A general 
solution of Eq. (17) is found by using the substitution 
7(t) 7(i) -I- k/Mrj in Eq. (2). In this case the variance 
of the process is equal to 



a2(t) = De-2fe*/Af7G2(i) 



GHs) 



ds. 



(18) 



The asymptotic stationary value of cr^(i) is = 
DMr]/2k — ksT/k which is independent from the pa- 
rameters a and f3. However, we observe a relaxation dy- 
namics whose functional form is controlled by the values 
of a and /3 parameters. To detect the different relax- 
ation dynamics we evaluate numerically the integral in 
Eq. (18) by setting 7(t) = a/{Tl^ + tl^). In Fig. 1 we show 
in a log- log plot the quantity A(t) = \<j'^ [t) — a1^\ / a% as a 
function of time. The quantity A(i) provides a measure 
of the distance of the system from the stationary behav- 
ior. In Fig. 1 we show that the quantity A{t) decreases 
following the power-law behavior A(i) cx for large 
values of t and for all the investigated values of the pa- 
rameters a and (3. In particular, when a > 0, <T^(t) — cr^f 
goes to zero as a negative value whereas when a < 
the same quantity goes to zero taking positive values. In 
order to illustrate this result we show cr^(t) as a func- 
tion of time when (3 = 0.6 and a = ±1 in the inset of 
Fig. 1. For < /3 < 1 (therefore including subdiffusive. 
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superdifFusive and stretched exponentially diffusing pro- 
cesses) it is not possible to define a characteristic time 
scale for the convergence of (j'^{t) during the process of 
relaxation. This is due to the fact that the integral of 
J^^ A{t)/A{ti)dt = oo. Although a power-law behavior 
is still observed when (3 > 1, it is worth pointing out 
that in this interval of P parameter a typical time scale 
might be determined by considering the above discussed 
integral which is finite in this region of /3 parameter. 

In conclusion the Langevin equations (1) and (17) with 
the choice of Eq. (8) describe non-stationary and station- 
ary random processes showing a wide class of (normal 
and anomalous) diffusion. When a stationary state ex- 
ists, the relaxation dynamics to the stationary state has 
a power-law time dependence. The processes modeled 
by Eqs (1) and (17) are characterized by a time depen- 
dent drift term in the associated Fokker-Planck equa- 
tion. Our model is complementary to the Batchelor's 
description of anomalous diffusion obtained by assuming 
a time-dependent diffusion term Equations (1) and^ 
(17) can be used to model metastable systems in which 
one of the physical observables, such as, for example, the 
viscosity, is time dependent. They can also be used to 
develop simple and efficient algorithms generating real- 
izations of random processes with controlled anomalous 
diffusion. 
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FIG. 1. Theoretical estimation of the normalized difference 
of the variance at time t from the stationary value a^, A(t) as 
a function of time. DiflFerent curves refer to different values of 
the control parameters a and j3. The a values are a = 1 (solid 
lines) and a = — 1 (dashed lines). The parameter (3 assumes 
the values 0.2, 0.6, 1.0, 1.4 and 1.8 from top to bottom. In the 
inset we show a typical time evolution of o-^(t) obtained by 
setting /3 = 0.6 and two different values of a, a = 1 (bottom 
curve) and a = — 1 (top curve). The other parameters are 
r = 1, k/Mr] = 1 and D = 1 

TABLE I. Summary of the different diffusion regimes. The 



constant C 


= 2|a|/(l 








a 


a\t) 


Description 


/3 


a = 


t 


Wiener (diffusive) 


/3> 1 
f3 > 1 


a > 
a < 


t 
t 


diflusive 
diflusive 


13^1 
(3^1 
f3 = l 


a > -1/2 
a = -1/2 
a < -1/2 


t 

tlogt 

^2|a| 


diffusive 
log divergent 
superdiffusive 


< /3 < 1 
< /3 < 1 


a > 

a < 




subdiffusive 
less than 
exponentially diffusive 


/3 = 
/3 = 


a > 
a < 


1 — exp(— 2at) 
exp(2|a|t) 


Ornstein-Uhlenbeck 
exponentially diffusive 


13 <0 
p <0 


a > 
a < 


1/tl^l 
exp[Cfi~'3] 


localized 
more than 
exponentially diffusive 
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